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Abstract 

^^- We study general features of three- meson bound states using the y(4660) as an 

(^ . example. Here the y(4660) is assumed to be either a two-body bound state of the 

O I /o(980), itself a bound state of K and K, and the ip' = V'(2'S), or a three-body bound 

state of ip' , K, and K. In particular, we investigate in detail the interplay of the 

various scales inherent in the problem, namely the /o binding energy, the Y binding 

energy, and the Ki/j' scattering length. This allows us to understand under which 

\^ ' circumstances the substructure of the /o(980) can be neglected in the description 

(N ! of they (4660). 
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O . 1. Introduction: In recent years a large number of new charmonium states 

has been discovered experimentally, most of them showing properties in vast 
conflict with the predictions of the quark model [1]. The structure of these 

/\ • states presents a serious challenge for our understanding of QCD. Since many 

c^ ■ of these new states are located near thresholds, it appeared natural to pro- 

pose a molecular nature for those. The prime example is the X(3872) located 
right at the D^D*^ threshold, proposed to be a molecule composed of these 
mesons [2] or at least a virtual state [3]. However, this picture is not uniformly 
accepted. It was, for example, challenged on the basis of data for X(3872) 
production in pp collisions at the Tevatron [1]. A rebuttal to this challenge 
was given in Ref. [S|. In the light quark sector, the light scalar mesons located 
below 1 GeV are also proposed to be molecules, namely the /o(600) — the a 
meson — the /o(980), ao(980), and k(900) [6|7] . Because of the relatively large 
widths involved, especially for the a, and the larger distance to the scatter- 
ing thresholds, however, the situation is less clear in this case. For a different 
viewpoint on the structure of these states, see 



These examples illustrate the importance of a model-independent method to 
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determine the structure of such states close to scattering thresholds. For two- 
body states located near a threshold there is a powerful method originally 
proposed by Weinberg for the deuteron [TU]. This method was extended in 
Ref. [TT] and allows one to quantify the two-body molecular component. A 
complementary formulation of this problem is given by effective range theory 
and the effective field theory (EFT) for large scattering length. The latter 
framework can also be extended to three- and four-body molecules [T2|T3] . 
Both methods effectively analyze the low-energy pole structure of the scatter- 
ing amplitude in a model-independent way. 

In this paper, we want to investigate to which extent it is possible to dis- 
tinguish between two- and three-body molecules and which scales govern this 
distinction. In particular, we analyze the example of the F (4660) in detail. This 
state was proposed to be a ip'f o{980) molecule based on an analysis with the 
Weinberg method [2] while the /o(980) itself was proposed to be a molecule 
of K and K mesons [SBj- It is therefore important to elucidate whether the 
internal structure of the /o(980) can be neglected for the description of the 
y(4660) or not. The answer to this question, of course, depends on the phys- 
ical scales in this problem. For a deeply-bound /o(980), one expects that the 
internal structure can be neglected. Performing numerically exact three-body 
calculations, we determine for which parameter ranges the 1^(4660) can be in- 
terpreted as a ip'f o{980) molecule and for which ranges as a ip'KK molecule. 
We stress that we only investigate the interplay of a possible two-body vs. 
three-body nature of the F (4660) in order to understand the interplay of the 
various scales. No explicit compact component is included for the F (4660) nor 
is the width of the /o(980) included. Thus, this work is mainly of theoretical 
interest and we are not yet in the position to compare with experimental line 
shapes. We start with a brief description of our formalism and some analyti- 
cal considerations before we present and discuss the results of our three-body 
calculation. The paper ends with an outlook. 

2. Formalism: In this section, we briefly outline the approach used to de- 
rive the amplitude for ?/^'/o(980)-scattering. For this purpose, we set up a 
non-relativistic effective field theory (EFT) for three distinguishable particles 
with different masseso Note that the typical momentum of the kaons in this 
problem is of order 70 MeV or less such that they can still be treated non- 



^ A detailed description of the EFT for three identical bosons is given in |13] . 
The general case of distinguishable particles with different masses can be found in 
Ref. [15]. 
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Fig. 1. Diagrammatic representation of the integral equations used. The bare and fuh 
dimer fields are indicated by the double and thick lines, respectively. The equation 
for the latter is illustrated in Fig. [2l 




Fig. 2. Diagrammatic representation for the full dimer propagator in channel k. 
relativistically. The underlying Lagrangian reads [15] 
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where the scalar functions ipk and dk are particle- and dimer fields, respectively, 
and rrik is the mass of particle k. The indices i,j,k in Eq. ([1]) are always 
different from each other such that the mass of particle k can also be labeled 
by rriij = rriji. This convention will be used below. Moreover, a two-body 
state of particles i and j can be labeled either by the indices ij or by the 
index k {i ^ j ^ k). The dimer fields are introduced for convenience and are 
not dynamical. An equivalent theory without dimer fields can be obtained by 
inserting the classical equation of motion for the non-dynamical dimer fields. 
Furthermore, Qk and h are coupling-constants, and the factors Ck sum up to 
one. From the Lagrangian ([T]), we can directly derive the Feynman-rules of 
our theory in momentum space and construct the general amplitude Aij for 
particle-dimer scattering in the channel i ^ j. The resulting equations are 
illustrated in Figs. [1] and [2j 

After renormalizing the amplitude and performing a partial wave decomposi- 
tion, the inhomogeneous integral equation for the /-th partial wave reads 
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where p and k is the incoming and outgoing momentum, respectively. The 



[Aj]i (P, k) = Y.r dg [T,k]i (p, q) Sk (g) [Ak,\ {q, k) - [T,,\ {p, k) , (2) 
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and the full dimer propagator 
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rrii + mj ' '" {rui + ruj) + ruk 

are the reduced masses of the two-particle system labeled by ij or by A; (z 7^ 
j 7^ k) and the particle-dimer system labeled by k, respectively. The Legendre 
functions of the second kind Qi and their arguments Cik{p, q) are defined by 
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where P; denotes the /-th Legendre polynomial. We use standard solution 
techniques in order to solve the integral Eq. (|2]) numerically. The a priori 
undetermined, dimensionless function H{A), that by construction contributes 
only in the s-wave, is directly related to the parameter h in Eq. ([T]). The func- 
tion H{A) plays the role of a running coupling constant and is determined by 
the renormalization procedure up to an unknown three-body parameter |16j . 
We fix its value by imposing that a three-particle bound-state exists at the 
binding energy E = —B. This requires that the amplitude Aij in Eq. ([2]) 
has a pole at this energy. However, any low-energy three-body observable can 
be chosen for this purpose. After fixing H, we are able to calculate the full 
particle-dimer scattering amplitude and therefrom, observables such as scat- 
tering lengths, cross sections and so forth can be predicted. The obtained 
results are independent of the value of the cutoff A, as long as A is large 
compared to all momentum scales in the problem. 

This theory represents the leading order of an EFT around the limit of large 
pair scattering lengths a^ [T21IT3] . Corrections are suppressed by powers of 
ro/ofc and pro, where tq is the range of the interaction and p is a typical mo- 
mentum. If the scattering lengths are not sufficiently large, the theory can still 
be taken as a particular three-body model for hadronic molecules. As such, 
it can be used to test the limits of applicability of the Weinberg approach to 
systems like the y(4660). By calculating the scattering length and the effec- 
tive range, we are able to investigate the question, whether the three-particle 
bound-state can effectively be treated as a two-particle system, consisting of 
one particle and an elementary dimer. If the substructure of the dimer be- 
comes relevant, the full three-particle picture is necessary. Transferred to our 



specific example y(4660), this means that we want to distinguish between the 
two alternatives ip'f q{980) and ip'KK. In the Weinberg formalism [10], the 
quantity Z G [0, 1] is defined as 
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where Y denotes the wave function of the physical state. Thus, Z measures 
the "non-two-particle" -fraction of the wave function. From the definition ([7]), 
we directly conclude that Z — >■ implies that the y(4660) is effectively a 
two-body system with an elementary /o(980)-dimer, whereas Z — )■ 1 means 
that it has to be seen as something else, in our case a three-body molecule. 
This method of distinction between two- and three-particle molecules can of 
course also be applied to other candidates for hadronic molecules in order 
to reveal their few-particle nature. The quantity Z is directly related to the 
residue Zpoie of the bound state pole in the two-particle scattering amplitude 
by Z = 1 - Zpoie- 

Performing a straightforward calculation within the Weinberg formalism, us- 
ing the Lippmann-Schwinger equation and the effective range expansion, the 
quantity Z can be related to the effective range parameters in two different 
ways [To]. The two corresponding Z-factors, which we denote by 

2P 

Za := °!? " :!:f° and Z. - ^^'^« 



V/o ~ ^'f'fo V/o 

should be (approximately) equal if the Weinberg formalism is applicable. The 
equations ([HI E]) receive corrections from the finite range of the interaction. 
They are exact in the limit of vanishing binding energy of the y (4660) relative 
to the -j/^'/o-threshold with Z = Za = Zr kept fixed [10]. The quantities a^'fo 
and r^//Q are the scattering length and the effective range for z/^'/o-scattering, 
respectively. Moreover, 

2P /'n^ 

'2 ^^if^B^if^ J 2 /i^,//Q {By - Bf^^ 



is the scattering length within a pure two-particle picture. 
The positive binding energies are defined via 

(10) 

(11) 
(12) 

Thereby all quantities on the right sides of the equations in ([8]) can be cal- 
culated in the full three-body model. Combining the two foregoing conditions 
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Z G [0, 1] and Z = Za = Z.,.^ we conclude that 

{)<Za^ Zr<l (13) 

should hold. Below we will use the validity of Eq. f ll3p as a diagnostic tool to 
identify the range of applicability of the Weinberg method. 

The full derivation of the equations ([H]) can be found in [10]. We omit it at 
this point, but outline its essential idea. Utilizing the binding energy between 
the particle ip' and /o(980) as defined in Eq. flTTl) . the Schrodinger equation 
can be applied in the form 

H\Y) = -B^,f,, I Y ) with H = Ho + V . (14) 

Here V describes the interaction and Hq denotes the free-particle part of the 
Hamiltonian, such that 

Ho\^'foip))=E{p)\^'fo{p)) (15) 

holds, where E{p) = p^ /{2^^if^) and p is the relative momentum of the if)' 
and /o in the center-of-mass. Employing the relations f lT^ and f lT^ . the 
momentum-integral appearing in the definition of Z (jTj) assumes the form 
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In order to deduce the formulas (jH]), the crucial requirement of the Wein- 
berg method is that the numerator on the right-hand-side of Eq. (fT6|) can be 
approximated by its value at zero momentum, 

, Mfo{p)\V\Y)\' ^ r \{^'MO)\V\Y)\^ 

For this to hold, the numerator has to vary slowly over the range of momenta 
contributing to the integral, i.e. 



Eip) < fi.0'/o ^ p^ ^^in^hB^o ■ (18) 

In other words, the range of the form factor {'ip' fo{p)\V\Y) has to be much 



larger than the range of the denominator, J2fi^rf^B^'f^y 

3. Analytical Considerations: Thus we have to understand the relevant 
momentum-scales of {'ip'foip)\y\Y). In the general case, with an incoming 
particle and dimer of type i this form-factor (i/'j di{p) \V\Y) is proportional to 
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Fig. 3. Feynman diagram for the form-factor {diilJi{p)\V\Y) with a particle and an 
dimer of type i in the incoming channel and y(4660) as the outgoing bound state. 

the bound state amplitude Fi as it is depicted in Fig. [3l It can be calculated 
by solving the homogeneous analog of Eq. ([2]) , 



[^^]o iP) = E /^ dg [T,,], (p, q) S, (q) [F,], (q) 
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where we have set / = for an s-wave state. In the rest frame of the F (4660), 
the total energy is simply E = —By- The quantity p is the relative momentum 
between the incoming particles. Solving Eq. (1191) exactly requires the same 
effort as solving the scattering equation and can only be done numerically. 
Such numerical solutions will be presented below. 

Our aim here is to investigate to what extent the underlying scales of the 
form factor {ipidi{p) \V\Y) can be understood in a simple way. For this pur- 
pose, we approximate the full amplitude [Fk]i (q) on the right-hand-side of 
Eq. ( TT9|) by a constant. In this point-like approximation, the p-dependence of 
{ipidi{p) \V\Y) emerges from the integral 
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where we have inserted the definitions ([3]) and (j4]). Using (Q, we deduce that 
for E = —By the quantity c := Cik{p,q) is always smaller than —1 and ap- 
proaches — oo for p — > 0. Hence a Taylor-expansion of Qq around 1/c = can 
be performed and we approximate: 
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Inserting this result into (l20l) . the momentum dependence of {ipidi{p) \V\Y) 
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so that the range of the form-factor can be estimated by A/2/ifci?y. Com- 
ing back to our specific situation with ip' and the /o-dimer, the range of 



W fo{p)\V\Y) is estimated as J2^^/kBy- From this we are able to formulate 
the condition 

yS/i^'E'-By ^ y 2/i^,//Qi?^//Q (24) 

for the validity of Eq. ( !T7|1 . Since the reduced mass terms /x^/.'^- and /i^'/o ^^^ 
both of the order rriK and -B,/,'/o = By — Bf^-^, this is equivalent to: 

BY:^BY-Bf^ ^ l>l-6o ^ 6o<l ^ --K1 ,(25) 

where we have introduced the dimensionless parameter Bq := Bf^jBy G (0, 1), 
which is the binding energy of the /o-dimer relative to the binding energy of 
the y(4660). Thus in this approximation the only relevant scale that controls 
the applicability of the Weinberg-method is bo which has to be close to 1, 
meaning that /o has to be a deeply bound state or an elementary particle. 

We stress that bo can not be the only relevant scale in this problem. In par- 
ticular, the scale characterizing the Kip' and Kip' interactions does not enter 
at all in this simple picture. In the full solution of the three-body problem 
where the y(4660) emerges as a dynamically generated three-particle-state, 
the situation will be more complex. We will investigate this question below by 
solving Eq. (|2]) numerically. 

4. Results and discussion: Using our formalism, we now present some re- 
sults for ■?/^'/o(980)-scattering within our three-body approach for the F (4660). 
In order to calculate observables, we first have to fix the 3 masses and scat- 
tering lengths, appearing in (E]), dl]) and (E]). Ignoring the - for our purpose - 
insignificant widths and errors, we take PDG values for the 3 masses by setting 
m^> = 3686.1 MeV and ttik = rrix = {fnK+ + 'tij^o)/2 = 495.6 MeV. The cor- 
responding two-particle scattering lengths are currently unknown, but for the 



kaon-kaon-interaction we can determine a^^ = ^1 y'^l^RRBfa by employing 
the analogue of Eq. ([9]). This corresponds to treating the /o as shallow bound 
state of K and K. For a^'K and a^//f, we have neither experimental data 
nor predictions from lattice calculations. Due to the absence of a two-particle 
bound-state and owing to symmetry reasons we can, however, conclude that 



a^p'K = a^'K < should hold. Another restriction on this quantity's magnitude 
conies from the fact that in '?/;'ii'-scattering no quark exchange is possible in 
contrast to D/C-scattering, for instance. Therefore, its interaction should be 
suppressed compared to the one in the latter case, implying that {a^'xl should 
be smaller than typical, non-resonant scattering lengths for Z^/T-processes. 
In [T7] , DK scattering lengths were calculated, using unitarised chiral pertur- 
bation theory and in Ref. [18] they were extracted from lattice QCD simula- 
tions. Both analyses agree and provide values of the order of 0.1 fm. Thus we 
demand —0.1 fm < a^'K = cl^'R < f^- 

We note that, given this assumption, a^'K is not large compared to the range of 
forces such that the EFT for large scattering lengths is strictly not applicable. 
In this study, however, we use the EFT as a specific three-body model that can 
be solved numerically in order to test the range of applicability of the Weinberg 
method. In relation to the corresponding thresholds in Eq. fITOl) and flT2l) the 
masses my = 4664 ± 11 ± 5 MeV and m^^ = 980 ± 10 MeV have large errors 
respectively so that we are still allowed to vary the binding energies in a region 
< Bfg < By ^ 20 MeV and be consistent with experiment. 




Fig. 4. Scattering length for the 'i/''/o(980)~process as a function of the relative 
binding energy bo = Bj^^/By for different values of By^ the binding energy of the 
Y with respect to the KKip' threshold (see Eq. (|11|) ). Our numerical values within 
the three-particle model (3P) are compared to the ones in a simple two-particle 



model (2P). The unknown scattering length a. 



ii'K 



ij^ is set to —0.1 fm. 



In Fig. m the numerically calculated '?/''/o(980)-scattering length within the 
three-particle model (3P) is depicted for several values of the total binding 
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Fig. 5. Effective range for '!/)'/o(980)-scattering as a function of the relative binding 
energy bo = Bf^/By for different values of By- The unknown scattering length 
Ciil)'K = clmk is set to —0.1 fm. 
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Fig. 6. Za from Eq. ([8]) as a function of the relative binding energy 6o = ^fo/^Y 
for different values of 6i = Bi/By. 
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energy By- It is also compared to the one of Eq. ([9]) in a hypothetical two- 
particle system (2P) with no substructure in the /o(980)-dimer. The variation 
in Bf^ is expressed in a dimensionless relative binding energy Bq. The un- 
known quantity a^>K = a^p'R is set to —0.1 fm. Our numerical results for the 
V''/o(980)-scattering length confirm the behavior deduced from the analytical 
considerations above (cf. Eq. fl25|l ). For a weakly bound /o(980)-dimer, that 
is bo — )■ 0, we see deviations between the two models, whereas for a strongly 
bound dimer with 6o — )■ 1 the results nearly coincide. In the latter case a^'f^ 
diverges to +00. We also calculated the effective range in Fig. Ofor the same 
parameters. It diverges at values where the scattering length vanishes. Un- 
fortunately, neither experimental data nor lattice calculations are available at 
present and no comparison to data for '?/;'/o(980)-scattering can be made. 

By applying the Weinberg-formalism, described in the previous section, we 
are also able to investigate the question at which scales F (4660) can effectively 
be seen as a two-body molecule consisting of elementary particles /o(980) and 
ip' and where the /o(980)-substructure has to be taken into account. Since 
the only undetermined scales in our system are By, Bf^ and a^'K, we expect 
the dimensionless quantities Za and Z^., as defined in the equations (jH]), only 
to depend on 2 dimensionless ratios of these 3 parameters. We have already 
shown above that the relative binding energy of the /o and y(4660), bo, plays 
an important role in determining the structure of the 1^(4660). For the second 
ratio it is natural to choose 61 := Bi/By := l/{2fi^iK cl^'k ^y)- Since the 
scattering length a^r^ is negative, the energy scale Bi characterizes a virtual 
state in the '?/''ii'-channel but not a bound state. The functions Za{bo) for 
discrete values of 61, varied over a wide range, are displayed in Fig. [61 We find 
that Za strongly depends on bo and 61. For 60 — ^ 0, Za reaches the value of 1 so 
that y(4660) has to be considered a three-body system. In the limit bo — > 1, 
Za approaches as the /o(980)-dimer gets maximally bound. In this case, the 
y(4660) can be seen as a two-particle system consisting of the elementary 
particles ip' and /o(980). 

Another feature shown in Fig. ([6]) is the fact that for certain parameters Za 
exceeds its allowed interval [0,1]. This leads us back to question at which 
scales the initial approximation (1171) is valid. The approximate condition for 
bo (1^ and the exact condition for bo and bi determined using Eq. (TT^ as 
a consistency-check for the theory, are summarized in Fig. [71 We show the 
consistency conditions for the Weinberg method in the {bo/bi)-{l/bo — 1)- 
plane. The shaded grey area indicates the allowed region with < Za < I 
while the area to the left of the dotted horizontal line satisfies Za ^ Zr- The 
area to the left of the dashed horizontal line satisfies the approximate condition 
(l25l) for 1/bo — 1^1. Within the lower left corner both our approximate 
and numerically exact calculations indicate the applicability of the Weinberg 
method. Outside of this area it leads to inconsistent results and can not be 
applied. For the case of the F (4660) this area is also mostly the region that is 
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Fig. 7. Consistency conditions for the Weinberg method in the 
{bo/hi)-{l/hQ — l)-plane. The shaded grey area indicates < Zq, < 1 while 
the area to the left of the dotted horizontal line satisfies Za ~ Zr- The area 
to the left of the dashed horizontal line satisfies the approximate condition for 
^0) ([25]) . Within the region of physical interest for the y(4660) (6i > 225), the 
consistency-conditions (|13p of the Weinberg-formalism are met almost anywhere. 



of physical interest. This is due to the fact that the experimental/theoretical 
constraints By ^ 20 MeV and |a^'x| ^ 0.1 fm restrict the parameter bi = 
l/{2fi^iK a^/K By) to be greater than approximately 225. 

5. Summary and Outlook: The work presented can not yet be used to 
compare to experimental line shapes, especially since we assumed the /o(980) 
as a stable particle. In reality, due the decay /o(980) — )■ vrvr the F(4660) 
becomes observable in i^'tttt invariant mass distributions and the resulting 
line shapes are believed to contain important information on the nature of 
the Y [13]. Within the formalism presented, this channel could be included 
via a complex scattering length of the KK system — the impact of unstable 
constituents on the line shapes of particles is, e.g., discussed in Ref. [TH]. In 
addition, the y(4660) may also decay into other channels. Currently there is 
a discussion, if the signal seen in A^A^, baptized X(4630), has its origin in 
the y(4660) [20]|21|22j — this channel could be parametrized by an imaginary 
part of the three-body interaction. 



However, we do not expect that the possible extensions mentioned will distort 
significantly the conclusions on the applicability of the Weinberg method. Our 
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analytical analysis summarized in Eq. ( 125|) shows that the binding energy of 
the y(4660) relative to the three-body threshold has to be large compared 
to the molecular binding energy relative to the ip'fo threshold. Moreover, we 
found that for systems where the scattering length in the other subsystems is 
large, this kind of analysis can not be used and should either be replaced by a 
more complex, coupled channel analysis or abandoned all together. However, 
in a significant part of the parameter space allowed for the /o(980) and the 
y (4660) the Weinberg method can be used to quantify the two-body molecular 
component, once better data is available. 

In principle, we could also integrate spin-effects, higher derivatives in the 
fields and, via photon-coupling, even charge-dependent interactions into the 
Lagrangian ([T]). Clearly all these extensions imply additional parameters in our 
theory which would have to be determined. Furthermore, we could scan the 
field of possible hadronic molecules for dimer- and three-particle-candidates. 
An interesting application is the X(3872)-meson as a DDn system. Since the 
Dvr-dimer could only appear in p-wave scattering, higher derivatives in the 
dimer fields would have to be included. 
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